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ABSTRACT 

Based on spectrophotometric observations from the Guillermo Haro Observatory 
(Cananea, Mexico), a study of the spectral properties of the complete sample of 24 
blue straggler stars (BSs) in the old Galactic open cluster M67 (NGC2682) is pre- 
sented. All spectra, calibrated using spectral standards, were re-calibrated by means 
of photometric magnitudes in the Beijing-Arizona- Taipei-Connecticut system, which 
includes fluxes in 11 bands covering ~ 3500 — 10000 A. The set of parameters was ob- 
tained using two complementary approaches that rely on a comparison of the spectra 
with (a) an empirical sample of stars with well-established spectral types and (b) a 
theoretical grid of optical spectra computed at both low and high resolution. The over- 
all results indicate that the BSs in M67 span a wide range in T c s (~ 5600 — 12600 K) 
and surface gravities that are fully compatible with those expected for main sequence 
objects (log g — 3.5 — 5.0 dex). 

Key words: stars: blue stragglers — stars: Hertzsprung-Russell (HR) diagram — 
Galaxy: open clusters and associations: individual: M67 



1 INTRODUCTION 

Blue straggler stars (BSs) were first discovered in the globu- 
lar cluster M3 by Sandage (1953). These peculiar stars were 
named 'blue stragglers' because of their observational prop- 
erties in star cluster colour-magnitude diagrams (CMDs). 
Usually, BSs appear as a bluer and brighter extension of 
a cluster's main sequence (MS). As members of the same 
star cluster and having been born at the same time, the 
behaviour of BSs is paradoxical because there should be no 
main sequence stars above the turn-off according to the stan- 
dard theoretical picture of stellar evolution in such a coeval 
and initially chemically homogeneous system. 

Decades have passed since their discovery, in which they 
have been the subject of many studies. These peculiar stars 
have been found to be common constituents of virtually all 
evolved systems (and also in young systems, but a 'normally 
populated main sequence' would hide any BSs), including 
dwarf galaxies (Stryker 1993). Based on observational and 
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theoretical studies, it is generally believed that the BSs in 
high-density regions of stellar systems could be the remnants 
of stellar collisions and those in sparse environments might 
result from the coalescence of interacting binaries or mass 
transfer through Roche-lobe overflow in primordial binary 
systems (Ahumada 1999; Bacon et al. 1996; Ferraro et al. 
1997; Gilliland & Brown 1992; Leonard 1989; Livio 1993; 
Ouellette & Pritchet 1998; Piotto et al. 1999; Stryker 1993; 
Tian et al. 2006). In addition to their still elusive origin, 
the study of BSs is important because in a stellar popula- 
tion they are among the most massive and luminous stars, 
whose contribution to the integrated light cannot be pre- 
dicted by the standard theory of stellar evolution (Bressan 
et al. 1993). In fact, it has been demonstrated that they 
greatly affect the spectral energy distribution (SED) of the 
entire population (Deng et al. 1999; Manteiga et al. 1989), 
particularly at ultraviolet and blue wavelengths (Xin et al. 
2007, 2008). 

In spite of the numerous studies published since their 
discovery, it is still not clear which of the conceivable expla- 
nations for the BS phenomenon is the preferred (or domi- 
nant) mechanism of formation. Similarly, it has not yet been 
established whether the spectral properties of BSs are the 
same as those of regular main sequence stars of the same 
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Table 1. The blue straggler population of M67. 'n' is the number 
of spectra collected for each object. 



Name R.A.(2000) Dec.(2000) ExpTime(s) n 

4 
2 
3 
4 
4 
2 
2 
4 
4 
2 
4 
2 
2 
2 
4 
2 
2 
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2 
2 
4 
4 
1 
2 



mass, as would be expected according to their loci in the 
CMDs, although this is in contrast to the potential chemi- 
cal enrichment in the atmospheres presumably provoked by 
the different detailed formation processes. 

Nevertheless, BSs have historically been regarded as 
core-hydrogen-burning stars (Benz & Hills 1987, 1992). For 
this reason, it is usually assumed that the spectral prop- 
erties of BSs are compatible with those of main sequence 
stars at the same loci in the CMDs. We have adopted this 
assumption throughout a recent series of papers discussing 
the integrated SEDs (ISEDs) of star clusters at low spec- 
tral resolution (Deng et al. 1999; Xin & Deng 2005; Xin et 
al. 2007; Xin et al. 2008). However, whether BSs can actu- 
ally be represented by main sequence objects has not yet 
been fully investigated, perhaps with the exception of rel- 
atively few early papers (Strom et al. 1971). The present 
paper is therefore aimed at validating this assumption ob- 
servationally. Determinations are also obtained of the two 
fundamental parameters, i.e., the effective temperature and 
the surface gravity, of the full sample of BSs in M67, based 
on a homogeneous collection of spectra. 

With the purpose of properly assessing their nature, 
we started a long-term project aimed at determining the 
atmospheric parameters of BSs in stellar systems. We will 
first determine the effective temperatures and surface grav- 
ities of the objects, through photometric and intermediate- 
resolution spectroscopic observations. In a second step we 
will investigate the chemical details for (corroboration of) 
a possible binary nature and to establish the existence (or 
absence) of the chemical patterns associated with a mass- 
transfer process. 

In this paper we present the initial steps of this project 
by investigating the full sample of BSs in the well-studied old 



Table 2. Artificial colours and magnitudes of our sample BSs. 
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open cluster M67 (NGC2682). M67 contains a rich system of 
24 BSs (Deng et al. 1999), a sample sufficiently large for sta- 
tistical purposes. The present paper is organised as follows. 
In Section 2 we describe the observations. In Section 3 we 
give the details of the flux-fitting method and provide the 
final sets of parameters. Fine-tuning of the gravity determi- 
nation is described in Section 4. In Section 5, a comparison 
with previous work is presented. Finally, a summary and the 
conclusions of this study are presented in Section 6. 



2 OBSERVATIONS AND REDUCTION 

The observations were carried out during a three-night run 
in February 2005 using the 2.12 m telescope of the Guillermo 
Haro Observatory (OAGH) at Cananea, Mexico. The spec- 
tra were collected using the Boiler & Chivens spectrograph 
with a 150 £/mm grating blazed at 5000 A and a Tektronix 
1024x1024 CCD detector. The instrumental set-up yielded 
a scale of 3.2 A per pixel with a wavelength coverage roughly 
from 3600 to 6900 A at a nominal 5.7 A full width at half- 
maximum (FWHM), with a slit width of 150 /jm. A total of 
66 object frames were observed, which included at least two 
frames per object, with the exception of BS206 for which we 
were able to observe only once. 

For the data reduction we followed standard procedures 
using the IRAF package. Bias and flat-field corrections were 
secured by collecting a set of ten bias frames as well as 
dome-projected halogen-lamp images at the beginning and 
end of each night. Each stellar image was accompanied by a 
Helium- Argon-lamp image that allowed wavelength calibra- 
tion and the determination of the nominal resolution along 
the dispersion axis. The relative flux calibration was done us- 
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Figure 1. Observed spectral energy distributions of the 24 BSs in our sample. The spectra are roughly ordered in a temperature sequence, 
decreasing from top to bottom. Vertical dotted lines indicate the loci of four major spectral features, He? 6562 A, H/3 4860 A, H7 4340 A 
and Ca K 3933 A. 
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Figure 2. Artificial colour-magnitude diagram for the BSs in 
M67. The full sample of 24 BSs are marked with open triangles. 
The solid line is the 4 Gyr isochrone of solar metallicity from 
Bertelli et al. (1994). 



Figure 3. A comparison of visual magnitudes (top panels) and 
B — V colour indices (bottom panels) of our programme stars 
between the present paper and those of Sandquist & Shetrone 
(2003) and Girard et al. (1989). 



ing observations of three spectrophotometric standard stars, 
BD75325, Feige67, and Feige34. 

The 24 BSs in our sample are all members of M67 with 
nearly 100 per cent membership probabilities, as determined 
from both proper-motion and radial-velocity observations 
(Girard et al. 1989; Sanders 1977). The catalogue is included 
in Table[T]where we give in columns (1) to (5) the BS identi- 
fication numbers from Fan et al. (1996), the equatorial coor- 
dinates, the integrated exposure times (in seconds), and the 
number of spectra collected for each object. The resulting 
relative-flux-calibrated spectra are shown in Fig. \T\ 

Qualitatively, the spectra in Fig. [1] are roughly ordered 
as a temperature sequence, based on visual inspection of 
the slope of the SED. It is interesting to note that the 
sequence of Balmer features, distinguishable down to Hll 
(A =3771 A), from bottom to top, exhibits an increase up 
to BS025, indicating that BS005 should have a temperature 
compatible with that expected for a late-B star. Similarly, 
the Ca K line at 3933 A is nearly absent in the two hottest 
objects and steadily increases in strength, overcoming the 
intensity of the blend with He and the Ca H line at 3968 A 
at the position of BS184. Another interesting feature easily 
observable in the spectra is the CH G-band at 4300 A. This 
feature is strongest for the two objects at the bottom. Again, 
from bottom to top, this feature disappears at the position 
of BS093. Therefore, this star should correspond roughly to 
spectral type A5 (about T cB = 8000 K). 

We apply an absolute-flux calibration using 
intermediate-band photometric data (resembling spec- 
trophotometric observations) collected at an earlier time 
(Deng et al. 1999). In brief, the 24 BSs were observed pho- 
tometrically using the Beijing- Arizona- Taipei-Connecticut 
(BATC) intermediate-band filters. These observations 
included 11 of the 15 filters in this system, covering a range 



from 3500 to 10000 A. The shortest wavelength covered by 
our dataset was obtained using a filter centered on 3890 
A. By convolving the observed BS spectra with the six 
intermediate-band (b,d, f, g,h, and i; central wavelengths 
at 3890 A, 4550 A, 5270 A, 5795 A, 6075 A, and 6660 A, 
respectively) transmission curves, six new magnitudes for 
each object were derived. These magnitudes were compared 
with those obtained with the BATC photometric observa- 
tions and permitted the derivation of the scaling factors 
to transform the spectroscopy-based magnitudes to the 
absolute BATC system. The accurate intermediate-band 
photometry in the BATC system secures (re-calibrates) the 
overall shape of the observed spectra for all programme 
stars. 

To assess the quality of the overall SED shapes and 
to check the precision of the calibration, the spectrophoto- 
metrically calibrated spectra were used to construct broad- 
band photometry, which can then be compared with stan- 
dard broad-band observations. Two sets of V, (B — V) obser- 
vations from Sandquist & Shetrone (2003) - who provided 
photometric data of BSs in M67 for a study of variability in 
the light curves - and Girard et al. (1989) - who studied the 
relative proper motions and the stellar velocity dispersion 
of M67 - are compared with the artificial magnitudes and 
colours (see Table [2| of the 24 BSs, derived by convolving 
the B and V filter-response functions with the calibrated 
spectra. Adopting a distance modulus of DM=9.97 nag and 
a colour excess E(B — V) of 0.059 mag for M67 (taken from 
WEBDAQ), the observations (both the location with respect 
to the cluster's main sequence turn-off and the magnitudes 
and colours of the programme stars) can very well be re- 
produced using the artificial CMD photometry, as shown in 

1 http://www.univie.ac.at/webda/ 
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Figure 4. Best fits for two representative BSs, BS018 and BS126. The solid and dashed lines show, respectively, the observed spectra and 
the best-fit Pickles (top panels) and Kurucz low- resolution spectra (bottom panels). The smaller panels below each spectrum correspond 
to the residuals, as explained in the text, and the dotted lines indicate the 3<r boundaries. 



Fig. [2] The artificial photometry was also compared with 
direct broad-band photometry (magnitudes and colours). A 
perfect match was found, as shown in Fig. [3] in which the 
residuals in V magnitude, 8 V=V-V t his work, and (B — V) 
colour index 5(B — V) = (B — V) — (B — V) t his work are dis- 
played. No systematic differences were found, whereas the 
random difference between the data derived from our spectra 
and from direct photometry is compatible with observational 
errors of a few per cent of a magnitude. Figure [3] indepen- 
dently shows that our spectral observations and calibration 
are accurate to a satisfactory degree. 



3 SPECTRAL FITTING AND ANALYSIS 

In order to study the spectral properties of BSs and to de- 
termine the effective temperature and surface gravity of our 



sample stars, we applied simple flux-fitting methods, using 
three different libraries of reference stellar spectra, both ob- 
served and synthetic. In all cases we assume a solar chemical 
composition for M67 stars, which is in agreement with ob- 
servational determinations (Bressan & Tautvaisiene 1996; 
Hobbs & Thorburn 1991). 

(i) Each calibrated spectrum was compared with every 
entry in the spectral atlas of Pickles (1998). The comparison 
was carried out after normalising our spectra and those in 
the reference atlas to the flux at A =5556 A. The algorithm 
we have implemented finds the spectrum (and its associated 
parameters) in the atlas that produces the minimum stan- 
dard deviation, a, of the residual flux, AF = Fbs — -Fpickies, 
computed for each A. 

Because of the marked decrease in sensitivity of the CCD 
at the shortest wavelengths, the spectral regime considered 
for the flux fitting excludes the region at A < 3850 A. Fig- 
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ure[4]displays the best fit and residuals for BS018 and BS126 
(top panels). The solid and dashed lines are, respectively, the 
BSs' calibrated spectra and the best-fit flux from the Pick- 
les library. The labels at the top indicate the star ID, its 
temperature and spectral-type designation. 

(ii) A similar procedure was followed by using the spec- 
tral grids of Lejeune et al. (1997, 1998) which are, for the 
segment of the parameter space under consideration, mostly 
based on Kurucz (1993) low-resolution theoretical fluxes. In 
this case, both sets of spectra (BSs and model fluxes) were 
normalised to the flux at A = 5390 A. A set of best-fit pa- 
rameters is found by directly comparing the observed spec- 
tra with each of the model fluxes. It is important to note 
that in this way, as well as in the previous point, the best 
fit always corresponds to a grid point. In Fig.|3]we show the 
best fit for the stars BS018 and BS126 [bottom panels). The 
solid and dashed lines are, respectively, the BSs' calibrated 
spectra and the best-fit theoretical flux. The label on the 
right gives the parameters of the best-fit model atmosphere. 

(iii) In this case we made use of the bluered library 
(Bertone et al. 2003, 2008). bluered is a high-resolution 
(i?=500,000) grid of over 800 synthetic stellar spectra, cov- 
ering SEDs in the optical range (A = 3500-7000 A). The 
library is based on the ATLAS9 model atmospheres and 
has been computed with the SYNTHE code developed by 
Kurucz (1993). The grid spans a large volume in the funda- 
mental parameter space, accounting for virtually any stellar 
type from O to M stars and from dwarfs to supergiants. 
An important aspect of this grid, although of marginal rel- 
evance for the parameters associated with our programme 
stars, is that its calculation includes the effect of diatomic 
molecules, in particular TiO. A best-fit spectrum was found 
in the two-dimensional space covering (T c s , log g) , after min- 
imising the statistical variance in the relative-flux domain 
as a measure of the similarity between target spectrum and 
theoretical SEDs across bluered (Bertone et al. 2004). As 
in the comparisons above, we have assumed a solar chemi- 
cal composition for M67. It is worth noting that the grid of 
theoretical spectra has been properly modified to simulate 
the instrumental set-up. The results are shown in columns 
(6) and (7) of Table [3] whereas contour plots for BS018 and 
BS126 are shown in Fig. [S] 

The results from the three methods are listed in Table [31 
where columns (1) to (8) include, respectively, the object ID, 
the parameter pairs (T e fr, log g or spectral type) and the 
identification numbers following Sanders (1977), for ease of 
cross identification. 

The agreement among the effective temperature esti- 
mates provided by the three methods is on the order of 2-6 
per cent and, in general, the best-fit T c ff values based on 
the bluered library are the highest, apart from the case of 
the hottest star, BS005, where the best fit is about 1550 K 
lower. The discrepancy in this case arises from the associ- 
ated low log g value, which is about 2 dex lower than the 
corresponding result from the Lejeune library, since the flux- 
fitting method, applied to intermediate-resolution spectra, is 
affected by a T c ff — log g degeneracy, where a lower surface 
gravity implies a cooler temperature (Buzzoni et al. 2001). 

A larger discrepancy affects the derived surface grav- 
ities of our sample stars, for which the highest values are 
most often provided by the bluered library. These latter 
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Figure 5. Contour plots for BS018 (top panel) and BS126 (bot- 
tom panel). The solid dot corresponds to the best parameter esti- 
mate; the contour levels indicate the 1, 2, 3, and 4cr uncertainties. 



results are, however, affected by an average la error of < 1 
dex (see Fig. [5]) The systematically higher parameter values 
of the bluered spectra can be understood as in Bertone et 
al. (2008), who show that the T e g and log g values that are 
obtained from comparing the Sun with the bluered spectra 
(at very high spectral resolution) are a few per cent higher 
than those commonly accepted because the physical param- 
eters of the absorption lines included in the spectral synthe- 
sis generate deeper features - which are counterbalanced by 
raising both the effective temperature and the surface grav- 
ity. In general, the current determinations of gravity for the 
programme stars are limited by the low resolution. 



4 FINE-TUNING OF THE GRAVITY 
DETERMINATION 

Complementary to the analysis presented in the previous 
section, we obtained observations at OAGH of the BS sample 
with an alternative set-up that allows, in principle, the sep- 
aration of potentially fiducial gravity indicators. In particu- 
lar, we will make use of the indicators defined by Rose (1984, 
1994), which consist of line ratios of several pairs of features 
and the corresponding index-index diagnostic diagrams, and 
the hydrogen- absorption indices defined by Worthey et al. 
(1994) as part of the Lick system. These two approaches are 
necessary in view of the large effective temperature interval 
covered by the BSs. 

New observations were carried out on February 24-27, 
2008, using the Boiler & Chivens spectrograph and the Ver- 
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Table 3. Best-fit parameters. 
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NOTE: a, Stellar Identification from Fan et al. (1996); b, Sanders number (Sanders 1977); *, Binary population. 

sarray 1300 x 1300 CCD detector optimised for the blue spec- 
tral interval. We used the 600 I /mm. grating and a slit width 
of 200 fim, which yielded a nominal dispersion of 0.7 A/pixel 
and a resolution of 2.6 A FWHM. The grating was posi- 
tioned to obtain spectra in the interval 3800-4700 A where 
all of the gravity indicators cited above are defined. 



The sample consisted of two stellar sets. The first corre- 
sponds to the full sample of BSs, whereas the other contains 
nearly 50 objects that served as gravity templates. The lat- 
ter set was selected from the catalogues of Cayrel de Strobel 
et al. (1997, 2001). Data reduction up to the standard flux 
calibration was performed using the conventional procedures 
of IRAF and utilising a set of standard stars observed each 
night. We considered it very important to secure calibrated 
fluxes to provide reproducible results when analysing data 
collected with other instruments. In Table U we list the con- 
trol sample, showing in columns (1) to (5) the stellar identi- 
fication, the spectral type and the associated stellar param- 
eters collected from Cayrel de Strobel et al. (1997, 2001). 
For the stars with multiple determinations we provide the 
average values. 



As an example of the flux-calibrated spectra we dis- 
play, in Fig. |6j the lower-resolution spectra of BS065 and a 
zoomed-in region at higher resolution. The vertical dashed 
lines indicate the position of several of the features used as 
gravity indicators, as described below. 




o.o t , , , u u , , 

3800 4000 4200 4400 4600 

wavelength (A) 



Figure 6. Low- (top panel) and intermediate-resolution (bottom 
panel) spectra of BS065. The dotted vertical lines indicate the 
positions of the Rose (1994) features used to define the gravity- 
sensitive line-depth index used in this paper. 



4.1 The wavelength sequence, line ratios and 
Lick-like indices 

We explored all possible combinations of the indices de- 
fined by Rose (1994), which are in the form of flux ratios 
at the central wavelengths of absorption lines or at pseudo- 
continuum loci, and the Lick/IDS indices (e.g., Worthey et 
al. 1994; Trager et al. 1998). We visually inspected the spec- 
tra of the template stars in search of additional pairs of 
features that could display a trend with gravity. At the end 
of the process, the indices that emerged as best gravity di- 
agnostics are (1) the combination 4289/4271 vs. H7/4325 
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Table 4. Stars used as gravity templates. 



Name Sp.Type T cS logg [Fe/H] 

(K) (dex) (dex) 

0.01 
0.40 
0.10 
0.14 
1.08 
0.30 
0.31 
0.23 
0.00 
0.16 
0.10 
0.07 
0.01 
0.06 
0.25 
0.22 
0.40 
0.17 
0.27 
1.46 
0.47 
0.04 
0.06 
0.22 
0.05 
0.78 
0.16 
1.00 
0.27 
0.25 
0.05 
0.38 
0.06 
0.05 
0.10 
0.05 
0.12 
0.52 
0.20 
0.18 
0.01 
0.20 
0.01 
0.05 
0.14 
0.48 
0.18 
0.23 
0.55 



from Rose's indices, for stars with T e g ^ 7500 K, and (2) 
the Lick index HyA, for stars hotter than 8000 K. 

4.2 Analysis of line-depth ratios 

As an important preliminary step, we theoretically verified 
the sensitivity of all of Rose's spectral features to surface 
gravity, in particular those that he identified as discrimina- 
tors of gravity. 

In spite of the similar spectral resolutions, this verifica- 




0.4 0.6 0.8 -5 5 10 15 

BLUERED index (A) BLUERED index (A) 



Figure 7. Comparisons between the theoretical and empirical 
indices for the sample of template stars. The empty circle in the 
right-hand panel corresponds to an object deviating more than 
3a. This object was not taken into account for the calibration. 
The dotted lines show the one-to-one correlations whereas the 
solid lines denote the best fits. 



Table 5. Linear transformation parameters. 



Index 


a 


b 


rms 


4289/4271 


0.385 


0.605 


0.017 


H7/4325 


0.120 


0.713 


0.026 


H 7A 


2.21 


0.80 


0.74 



tion process is needed because Rose's spectra were not flux 
calibrated and, therefore, subject to effects inherent to the 
particular instrumental set-up he used. In other words, we 
have corroborated that the selected indices indeed separate 
the effects of effective temperature from those of gravity. 

We calculated the line-depth ratios in a subsample of 
solar chemical composition synthetic spectra of bluered, 
after properly degrading the grid to match the working res- 
olution of 2.6 A. As mentioned previously, after exploring 
the full set of line ratios, we identified the diagnostic dia- 
gram including 4289/4271 vs. H7/4325 as the best fiducial 
combination for differentiating amongst stellar luminosity 
classes. 

Once we have confirmed their sensitivity to gravity, the 
next step is to transform the theoretical indices to our ob- 
servational system. For this calibration we compare the the- 
oretical indices with those measured from observed spectra. 
Note that theoretical values were obtained from bluered 
spectra after linearly interpolating the set of parameters of 
the template stars (Table 2} . As an example of this com- 
parison we show in Fig. [7] the correlations between the the- 
oretical and the empirical indices for R~7/4325 and H7A. 
This figure indicates that a linear transformation of the form 
indexthaoi = a + bx index h s , suffices to properly match the 
theoretical and empirical indices. 

The comparison resulted in the transformation coeffi- 
cients listed in Table [5j along with the root mean square 
error. 

In Fig. [S] we display the theoretical diagram (for solar 
metallicity) for a set of different effective temperatures and 
gravities, after application of the transformation described 
above. The dashed and solid lines illustrate, respectively, the 
iso-Tgff and iso-gravity curves. In the top panel we overplot 
the loci of the template stars with —0.15 ^ [Fe/H] <+0.15 
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Figure 8. Diagnostic diagram for Rose's indices 4289/4271 vs. 
H 7 /4325. Dashed lines indicate the iso-T G ff curves at 6000, 6500, 
7000, 7500, 8000, and 8500 K (from right to left), while the solid 
lines show iso-gravity trends at logg=5, 4, 3, 2 and 1 dex for 
the theoretically calculated indices. In the top panel the upward 
and downward triangles indicate, respectively, reference stars with 
gravities in the intervals log g=4. 0—4.5 dex and log g=3. 5-4.0 dex. 
In the bottom panel the positions of the BSs are marked with open 
circles, along with their identifications. 



dex. In the bottom panel we show the same diagram and 
the positions of the BSs. 

From inspection of the panels in Fig. [8] we note the 
following characteristics: up to an effective temperature of 
7500 K, the indices can clearly be used to separate among 
stars of the three higher gravity bins, a sensitivity which ap- 



pears enhanced for dwarfs and subgiants. Most of the BSs 
show up in the interval log g=4. 0-5.0. The only exceptions 
are BS034, BS043, BS065 and BS216, for which their loci in 
the diagrams indicate gravities lower than logg = 4.0 dex. In- 
terestingly, in the work of Mathys (1991) two of these stars, 
BS034 and BS043, also turned out to be the lowest-gravity 
objects, with logg=3.79 and 3.44, respectively, although the 
latter star might require a more detailed analysis since it is 
part of a close pair (Girard et al. 1989). 

Importantly, we note at this point that we have not 
(yet) attempted to use these diagrams to derive values for 
the atmospheric parameters, but instead we only provide an 
overall assessment of the gravity of the objects. For a more 
quantitative evaluation, a detailed analysis regarding the ad- 
equacy of the model spectra is necessary and still beyond the 
scope of the present paper. At any rate, it is important to 
mention that in spite of the potential problems associated 
with theoretical spectra (see Bertone et al. 2008) the dia- 
grams clearly exhibit (even without prior knowledge of the 
atmospheric parameters) that most stars with T e g ^ 7500 K 
have surface gravities compatible with stars on the main se- 
quence. 



4.3 The H/y Lick-like index 

For the BS stars of higher temperatures we have measured 
absorption indices of the hydrogen Balmer lines, as denned 
by Trager et al. (1998). We have termed these indices 'Lick- 
like' since they have not been transformed to the Lick sys- 
tem. The overall behaviour of the indices associated with 
the H7 and H<5 lines in empirical data has demonstrated 
that the indices barely depend on metallicity and are very 
sensitive to gravity for stars with Terr > 8500 K. We have 
constructed a theoretical diagnostic diagram of H7A vs. T e g 
using bluered. In a similar fashion to the line-depth in- 
dices, we have calibrated theoretical indices by comparing 
them to the indices measured in the template stars. The 
results are included in Table(5] In Fig. |9]we display the the- 
oretical trends for solar chemical composition, with the stars 
represented as in Fig. [5] In Fig. |5]we include the full sample 
of template stars regardless of their chemical composition. 
Note that the index values degenerate at low temperatures, 
whereas stars are well separated at high temperatures, in 
particular BS005 (our hottest object). According to this di- 
agram, the three hottest stars have surface gravities in ex- 
cess of log g=4.0, although - because of their temperature - 
their loci in the diagram do not allow us to precisely estab- 
lish this parameter. The hot object BS005 appears to have 
a gravity of about log<j=3.6, which is compatible with our 
determination using the bluered grid. 

There are two objects, BS029 and BS038, that do not lie 
within the physically expected regions in the two diagrams. 
For these stars, Mathys (1991) determined effective temper- 
atures consistent with our determination, and gravities of 
log 3=3.91 and 4.14 for BS029 and BS038, respectively. 

Therefore, the spectral properties of BSs can be repre- 
sented by main sequence stars with the same photometric 
properties when modeling a simple stellar population based 
on (photometric) observations of star clusters. 
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Figure 9. Lick-like index H7A as a function of effective temper- 
ature. The solid lines show the theoretical iso-gravity curves from 
log g=5 to 1 dex from solar metallicity spectra. Observed stars 
are marked with the same symbols as in Fig. [8] with the addition 
of squares marking objects with log 3=2.5—3.5 dex and starred 
symbols denoting stars with log g <2.5 dex. 



5 COMPARISON WITH PREVIOUS WORK 

Previous studies of BSs in M67 include, for instance, Bruntt 
et al. (2007) for BS018, BS025, BS034, BS038, BS040, 
BS047, and BS093, based on asteroseismic analysis for 
5 Scuti pulsations, and Shetrone & Sandquist (2000) for 
BS043, BS046, BS139, and BS206, based on abundance anal- 
ysis. Mathys (1991) analysed 11 BSs in M67 and studies on 
binarity are also available (see below). However, a homoge- 
neous survey of the full sample of BSs in M67 and complete 
atmospheric parameter determinations have not yet been 
carried out. Mathys (1991) presented a spectroscopic study 
of 11 BSs in M67, and performed a detailed abundance anal- 
ysis for F 153 and F 185. He concluded that the effective 
temperatures and surface gravities of the BSs in M67 were 
quite similar to those of normal main sequence stars of the 
same spectral type. 

There is an obvious difference between his method 
and ours. Based on the photometric data from Mermil- 
liod (1988), Mathys (1991) derived the atmospheric param- 
eters from the photometric measurements of the BSs in the 
Stromgren system, applying the relevant calibration to de- 
termine the effective temperatures and surface gravities of 
B-, A- and F-type stars using uvbyfl photometry (Moon & 
Dworetsky 1985). 

The effective temperatures derived from the present 
work (from bluered) for the 11 BSs in common are con- 
sistent with Mathys (1991), as shown in Fig. HOf a^l . and 
the surface-gravity determinations are less conclusive for 
most of the BSs compared with Mathys (1991), as shown 
in Fig. HOf bh The error bars (horizontal axes) in Fig. 1101 
were obtained based on bluered spectral fits, whereas the 
vertical error bars are from Mathys (1991). 



For the 11 BSs in common, the surface gravities of 
BS005, BS018, BS025, BS038, BS040, BS046, BS047, and 
BS093 in both papers are in fairly good agreement, consid- 
ering the error bars. Large deviations in surface gravities 
are found for the remaining three BSs in common (BS029, 
BS034, BS043), as shown in Fig-E^b). Very likely, the un- 
doubted binary nature of these three objects is responsible 
for the deviations between the two methods. 

Indeed, there are eight objects in the list of BSs in M67 
which are likely binary candidates, based on previous ob- 
servations. The BSs identified as binaries are marked by as- 
terisks in Table H The BS BS034 (S1284), for instance, is 
thought to be a binary system in the final stages of mass 
transfer (Milone & Latham 1992; Zhang et al. 2005). Milone 
& Latham (1992) considered the dominant light contributor 
in BS034 to be the original primary (now the secondary) 
with an orbital period of 4.18284 days and an eccentricity 
of e=0.205. Based on high-precision radial velocity measure- 
ments, they obtained a spectroscopic orbital solution for the 
BS binary system. They supposed that the mass transfer be- 
gan fairly recently and that this BS was formed through sta- 
ble mass transfer with nearly 100 per cent efficiency. BS046 
(S1082) was found to be a complex unusual eclipsing binary 
system, or even a triple system of which the SED could be 
explained by the sum of a close binary and another main 
sequence star (van den Berg et al. 2001; Zhang et al. 2005). 
BS029 (S1267), BS043 (S975), BS047 (S752), BS111 (S997), 
and BS115 (SI 195) were all identified as spectroscopic bi- 
naries with long periods ranging from 800 to 5000 days 
(Latham & Milone 1996). BS184 (S1036) was detected as 
a W UMa-type binary with a small amplitude of light vari- 
ations, and a strong but stable O'Connell effect (Sandquist 
& Shetrone 2003). 

In our work, these binaries were easily fitted using 
model spectra of single stars. Although we cannot corrobo- 
rate their binary nature, to within the limited resolution of 
the observations, we can nevertheless provide constraints on 
the BSs' spectral properties. 



6 SUMMARY AND DISCUSSION 

This study represents the first attempt to derive the param- 
eters of the full sample of BSs in the old Galactic open clus- 
ter M67 (NGC2682) in a homogeneous way. Low-resolution 
spectra of the sample of 24 BSs in M67 were collected us- 
ing the 2.12 m telescope of the Guillermo Haro Observatory 
(Mexico). The entire data set was re-calibrated using the 
BATC intermediate-band photometric system, in addition 
to the usual relative calibration using standard stars, and 
was subsequently used for a comparison with three different 
stellar databases aimed at studying their spectral properties 
in a systematic way. We found that all objects have gravity 
values in agreement with the expected values for objects in 
the hydrogen-burning stage. 

Considering the original goal of our work, we conclude 
that, in terms of spectroscopic properties at low resolution, 
the BSs can indeed be represented by empirical or theoret- 
ical data of (or compatible with) main sequence stars, at 
least in a low density environment as in M67. 

As a natural extension to this, it is further concluded 
that when building up the empirical SEDs of SSPs based 
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Figure 10. Comparison between our results and those of Mathys 
(1991): (a) effective temperature; (b) surface gravity. The dotted 
lines are the one-to-one correlations. The open circles are the 11 
BSs (labeled by star IDs) in common between the present work 
and Mathys (1991). 



on stellar clusters, the contributions due to BSs can be ac- 
counted for using photometry and stellar spectral libraries. 
This conclusion holds at least at low and intermediate spec- 
tral resolution. 

Limited by the spectral resolution of the current obser- 
vational data set, it is not possible to assess binarity and 
the formation mechanism of the sample of BSs in M67. We 
anticipate that a detailed chemical abundance analysis at 
high resolution will show signatures of these dynamical and 
physical processes. Therefore, the current work serves as a 
valuable starting point. 
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